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Oscillatory behavior and enhancement of the surface plasmon 
^ : linewidth in embedded noble metal nanoparticles 
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Abstract. We study the Landau damping of the surface plasmon resonance of metallic nanoparticles em- 
bedded in different environments of experimental relevance. Important oscillations of the plasmon linewidth 
as a function of the radius of the nanoparticles are obtained from numerical calculations based on the time 
dependent local density approximation. These size-oscillations are understood, within a semiclassical ap- 
proximation, as a consequence of correlations in the spectral density of the nanoparticles. We treat inert 
matrices, as well as the case with an unoccupied conduction band. In the latter case, the plasmon lifetime 
is greatly reduced with respect to the inert case, but the non-monotonous size-dependence persists. 
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1 Introduction 

In recent years, the use of femtosecond laser pulses has 
made it possible to investigate the time evolution of col- 
lective and quasi-particle excitations in metallic nanopar- 
ticles. In these experiments, the actual time-resolved evo- 
lution of the energy transfer to the environment can be 
studied A determining factor in the relaxation pro- 

cesses is the lifetime of the collective excitations (surface 
plasmons), which for small cluster sizes (0.5 nm < R < 

2 nm) , is limited by Landau damping (decay into electron- 
hole pairs). The lifetime of the surface plasmon resonance 
is thus, not only a problem of fundamental interest, but 
also important for possible applications. 

In the last decade, the analogy of surface plasmons 
with nuclear giant resonances has been exploited, in order 
to obtain a physical understanding, as well as analytical 
and numerical estimations, of the plasmon lifetime 
0. Following these approaches, and using a semiclassical 
formalism, an oscillatory behavior of the level width was 
determined for free clusters of small sizes (a number of 
atoms N between N = 20 and N = 1000). These analyt- 
ical results have been verified against numerical calcula- 
tions within the time-dependent local density (TDLDA) 
for a jellium model ||], and the oscillations in the level 
width with the size of the nanoparticles were shown to 
arise from the electron-hole density-density correlations 
in the angular restricted density of states. 

The non-monotonous size-dependence of the plasmon 
lifetime has been experimentally observed in free alkaline 
metal nanoparticles j!|l(J, as well as in embedded noble 
metal clusters [fllj. Further experimental and theoretical 



work would be needed in order to unambiguously char- 
acterize the level width oscillations. While the analytical 
calculations are most easily done for free alkaline clusters, 
metallic noble metal clusters embedded in a transparent 
matrix are more conveniently investigated from the exper- 
imental point of view Q . It is then important to study the- 
oretically the size-dependence of plasmon lifetime for em- 
bedded particles. This is the purpose of this work, where 
we extend the results of Ref . [ 7| to different kinds of envi- 
ronments. 

The problem of the plasmon lifetime in metallic nano- 
particles embedded in a matrix, is considerably more diffi- 
cult, and less understood, than in the case of free clusters. 
The most obvious effect when the clusters are embedded 
in a dielectric material, is the change produced by the di- 
electric constant of the matrix, that moves the position of 
the plasmon resonance according to the Mie formula jl3| , 
and can also affect the width of the plasmon jjj. Moreover, 
when the matrix is not inert, it can react with the atoms 
in the surface of the nanoparticle. The electronic states 
on the surface are modified and contribute to the loss 
of coherence of the collective state. This is the so-called 
chemical interface damping. In addition, the presence of a 
conduction band in the matrix at an energy accessible to 
the electrons, weakens the confinement, and contributes to 
the width of the surface plasmon in an important way. All 
these different effects can add up and make the theoretical 
description of the system quite involved. 

The experiments of Charle et al. |llj studied the in- 
fluence of the matrix environments on surface plasmon 
excitations of small silver particles. They found a strong 
broadening of the surface plasmon resonance when the 
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particles were embedded in a reactive matrix such as CO, 
as compared to an inert matrix from a noble gas like 
Ar. This effect has been attributed to chemical interface 
damping. Also, a very strong broadening has been ob- 
served by Hovel et al. ]l4j for silver particles embedded 
in SiC>2. In a very interesting and pioneering paper, Pers- 
son proposed a model for explaining these experimental 
observations that took into account the modification of 
the electronic eigenstates of the metallic nanoparticle by 
the layer of absorbates on its surface and by a conduction 
band |L§. 

In the next section we recall the basic theoretical facts 
of the calculation of the plasmon lifetime within a semi- 
classical approximation. We then study the linewidth of 
embedded clusters from TDLDA calculations for inert ma- 
trices. The case of matrices with a conduction band is 
treated in Sec. ||. The theoretical treatment of the chem- 
ical interface damping is out of reach of our theoretical 
description, but the cases that we will consider are not 
thought to be considerably influenced by this effect. 



2 Theoretical background for free clusters 

The Landau damping of the dipole plasmon can be calcu- 
lated treating the collective excitation as an external per- 
turbation, which can give rise to the creation of electron- 
hole excitations |l2| . Then, Fermi's golden rule yields the 
linewidth 



Ar(R) = 2nJ2\(p\SV\h)\ 

ph 



5(hid M - e P + e/i) 



(1) 



where \p) and \h) are electron and hole states in the self- 
consistent field, with its energies given by e p and e/j, re- 
spectively, 5V is the dipole field due to the surface plas- 
mon, and wm is the frequency of the plasmon, which clas- 
sically is given by the Mie formula, 



\feJT2e~r. 



(2) 



where w v is the plasma frequency, td is the dielectric func- 
tion of the d electrons and e m is the dielectric function of 
the embedding matrix 

In the case of spherical symmetry, assuming that the 
confinement and the interactions lead to hard walls at 
radius R in the self-consistent field, we can evaluate Eq. 
CPT) . Integrating over the electron- hole states, one obtains 



AT(R)=cry f dEj2 E (2L + 1)(2L' + 1) (3) 



L L' = L±1 



x < L, 0; 1, 0|£', > 2 E(E - haj M )d L (E)d L ,(E - Tiw 



M J 



where < L,0; 1,0|L',0 > is a Clebsch-Gordan coefficient, 
7 = (2irh 3 )/{3NM 2 u) M R 4 ), c = 4Mi? 2 /S 2 , E F is the 
Fermi energy, and dh{E) is the one-dimensional density 
of states with total angular momentum L. 



Using the semiclassical expression for the density of 
states of the one-dimensional problem |l6[| , we can decom- 
pose d^E) in its smooth (zero-length trajectories) and 
oscillating components (arising from periodic orbits). We 
then obtain two contributions to the width of the plasmon 
resonance R 



AF = r Q + r a 



(4) 



The smooth term Fq arises from the smooth component 
of d^ and exhibits the well known 1/R dependence firstly 
proposed by Kawabata and Kubo(l7|,|j| 



(5) 



where £ = Hujm/Ey , <?(£) is a smoothly decreasing func- 
tion with g(0) = 1, and vf is the Fermi velocity. 

The oscillating part r osc arises from the density oscil- 
lations as a function of the energy. Within a semiclassical 
approach, it can be written as R 



6V27T>1 



MR 2 y/k^R^ ~J \pr 



E ~r cos(2rfc Fj RO, (6) 



where the sum runs over all repetitions r of the period of 
the equivalent one-dimensional motion. The amplitude of 
this oscillations can be of the order of r for small clusters. 



3 Embedded metallic clusters in inert 
matrices 

We want to address the question of what happens with the 
oscillations and the typical value of the linewidth of the 
plasmon when the nanoparticles are embedded in a ma- 
trix. In Ref. § we presented some calculations for noble 
metals embedded in inert matrices. The numerical results 
were obtained using the TDLDA formulation by Bertsch 
p8| but modifying the residual interaction to include the 
frequency-dependent dielectric function of the d electrons 
ed(w) and the dielectric constant of the matrix e rn [ p"9| . In 
order to maintain the spherical symmetry of our problem 
we always considered cluster sizes corresponding to magic 
numbers of atoms. 

For inert matrices, the plasmon linewidth is affected 
only by the modification of the frequency lom through the 
parameter £ of Eqs. (||) and (||) . This is due to the fact that 
the self-consistent potential, and, therefore, the electron- 
hole density-density correlation function are practically 
unchanged by e m , but the energy position of the plas- 
mon peak follows the Mie formula. Results for the width 
of the plasmon in the case of Ag nanoparticles embedded 
in a matrix of Ar are shown in Fig. [l| In the inset we show 
a typical spectrum of the photo-absorption cross-section 
for Ag44o embedded in Ar. The singularities of the spec- 
trum are smeared out by a non-zero 7. This value of 7 
is subtracted at the end of the calculation and we verify 
that the fit of AF is not sensitive to it, provided that it 
is larger than a minimum value. Ar exhibits pronounced 
oscillations as a function of the radius. The smooth part of 
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Fig. 1. Linewidth as a function of the radius (in units of 
the Bohr radius ob = 0.53 A) for silver nanoparticles in an Ar 
matrix calculated within TDLDA (full circles) together with 
the experimental results of Ref. [[HJ (empty squares) . The dot- 
ted line through the numerical points is a guide-to-the-eye. 
The dashed and the solid lines represent, respectively, Jo(i?) 
and Ar(R) according to Eqs. (^), (|B|) and (fl) (with a reduc- 
tion factor of 3 as discussed in the text). Inset: logarithm of 
the TDLDA absorption cross section for Ag44o (R/as = 2.2), 
showing the pronounced surface plasmon resonance, fitted by 
a Lorentzian (dotted line). The excited states are indicated by 
tick marks and their oscillation strengths given by the height 
of the vertical lines. 



the linewidth, and our semiclassical result, are an overall 
factor of 3 larger than the numerical results. Therefore, 
they have been rescaled for in the figure. The experimen- 
tal values of Charle et al. (empty squares) are relatively 
well described by the TDLDA results (within a 20 %). The 
oscillations are suppressed due to the large particle sizes 
and the smearing resulting from wide size-distributions. 
The difference between the Kubo formula, on one hand, 
and the numerical and experimental results, on the other 
hand, is not understood. In any case, given the approx- 
imations used in the semiclassical calculation, it is not 
surprising that, although, this simple theory works well 
with alkaline metals, when we include e^w) there are fac- 
tors not taken into account. Nevertheless, the analytically 
obtained oscillations have the correct period and relative 
size. 

Moving through different values of e m we change the 
position of the plasmon frequency ujm- Since the electron- 
hole density-density correlation function is insensitive to 
e m , the plasmon lifetime behaves in a non- monotonous 
fashion as a function of e m . Thus, we can use the dielec- 
tric constant of the embedding medium as a probe for 
this density-density correlation function. As an example, 
we show in Fig. || the results of TDLDA calculations for 
the width of the surface plasmon resonance for silver clus- 
ters of two different sizes (N — 92 and N = 198) as a 
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Fig. 2. Width of the surface plasmon resonance for silver clus- 
ters of two different sizes as a function of e m . Inset: frequency 
of the plasmon ujm together with the expression for the Mie 
formula (solid line) with a small modification of the bulk ui p . 



function of e m . The oscillations in the width are very dif- 
ferent for the two sizes due to the shell filling and the 
corresponding density correlation function. Although the 
range of e m over which we can see the complete oscilla- 
tion of the plasmon lifetime is very large, it should be 
possible to observe maxima and minima in experiments, 
if we are able to change the dielectric constant of the sur- 
rounding medium. In the inset of this figure we show the 
position of the plasmon peak as a function of the e m for 
the corresponding sizes. The behavior with e m is very well 
described by the Mie formula, and the two sizes behave in 
the same way. We can see almost no difference in the po- 
sition of the plasmon peak between the sizes for the same 
value of e m , but there are great changes in its width. 



4 Ag nanoparticles embedded in a matrix 
with a conduction band 

The SiC>2 used in experiments for embedding Ag nanopar- 
ticles is an amorphous solid with a conduction band with a 
minimum situated at — 1.7 eV with respect to the vacuum 
energy. The valence band maximum occurs 10.6 eV below 
the vacuum energy and has no influence in the Ag-surface 
plasmon |l5| . Chemical interaction in the surface between 
Ag and the SiC-2 is not expected to occur iQ and chemical 
interface damping does not influence the width of the plas- 
mon in this case. In order to implement the TDLDA cal- 
culations, we simulate the embedded medium by a change 
in the boundary conditions for the calculation in the self- 
consistent potential V(r) so that V(r) — > — 1.7eV when 
r — ► oo. In this way the electrons are less bounded to 
the cluster and can go to the conduction band in the ma- 
trix. The results of the self-consistent calculation of V(r) 
and its comparison with the potential without conduction 
band are shown in the inset of Fig. 0. We can clearly see 
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that the electrons at Ep are less tight, which translates in 
a small redshift in the position of the plasmon peak due 
to the bigger spill-out and into an increase (by a factor of 
2) of the width (compare with Fig. |l|). This is consistent 
with the factor of 3 enhancement observed in the experi- 
ments of Ref. Il4|] . Although in these experiments the dis- 
tribution of sizes has certain unknown variance that can 
increase the results for the linewidth of the optical ab- 
sorption experiments, a more refined model including the 
self-interaction correction to the TDLDA, and a corrected 
dielectric function for the surface, are probably needed for 
a more quantitative agreement. However, we can clearly 
see that the non-monotonous behavior is maintained and 
is of the same order as for free Ag nanoparticles. Other dif- 
ferent kinds of glasses that are used in experiments should 
produce similar increases in the width of the plasmon if 
the difference between the energy of the minimum of the 
conduction band and Ep is not too far from the energy of 
the plasmon. 



5 Conclusions 

As other size-dependent phenomena in clusters, the width 
of the plasmon presents a leading-order (smooth) contri- 
bution (that goes like 1/R), with oscillatory corrections 
due to shell effects. For the plasmon lifetime these size- 
dependent corrections are much more important than for 
the position of the resonance. Going for free to embedded 
clusters makes this difference still more pronounced. 

We have shown in this work that the size-dependent 
oscillations of the linewidth, in free as well as in embed- 
ded clusters, arise from electron-hole density correlations. 
Such an effect can be explored by varying the dielectric 
constant of the matrix around the nanoparticle. 

We have also considered the width of the plasmon for 
a simple model of nanoparticles of Ag surrounded by a 
matrix of amorphous SiC>2 with a conduction band. The 
width is increased in a way that agrees semi-quantitatively 
with the experiments. The oscillations in the width are 
still present and could be seen in experiments with a nar- 
row distribution of sizes. An adequate choice of the com- 
position of the matrix should make it possible to increase 
the lifetime of the surface plasmon by tuning its energy 
to a value where the electron-hole correlation function is 
nearly zero, which can be useful in certain applications. 

We thank J.-Y. Bigot, H. Haberland, and P.A. 
Hcrvicux for helpful discussions. RAM acknowledges the 
financial support from the European Union's Human Po- 
tential Program (contract HPRN-CT-2000-00144). 
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Fig. 3. Width of the surface plasmon resonance as a function 
of the radius of the particle for silver clusters embedded in SiC>2 
(the matrix is assumed to support a conduction band, see text). 
The discontinuous line is the best fit to 1/72 behavior. Inset: 
comparison of the self-consistent potential V(r) with (full line) 
and without (dashed line) the conduction band. We indicate 
the position of the Fermi energy with an arrow. 



3. T. Klar, M. Perner, S. Grosse, G. von Plessen, W. Spirkl, 
J. Feldmann, Phys. Rev. Lett. 80, 4249 (1998). 

4. M. Barma, V. Subrahmanyam, J. Phys.: Condens. Matter 
1, 7681 (1989). 

5. C. Yannouleas and R.A. Broglia, Ann. Phys. (N.Y.) 217, 
105 (1992). 

6. C. Yannouleas, E. Vigezzi, R.A. Broglia, Phys. Rev. B 47, 
9849 (1993); C. Yannouleas, Phys. Rev. B 58, 6748 (1998). 

7. R.A. Molina, D. Weinmann and R.A. Jalabert, Phys. Rev. 
B 65, 155427 (2002). 

8. W. Ekardt, Phys. Rev. B 31, 6360 (1985). 

9. C. Brechignac, Ph. Cahuzac, J. Leygnier, A. Sarfati, Z. 
Phys. D 12, 2036 (1993). 

10. S. Mochizuki, M. Sasaki, R. Ruppin, J. Phys.: Condens. 
Matter 9, 5801 (1997). 

11. K.P. Charle, W. Schulze, B. Winter, Z. Phys. D 12, 471 
(1989). 

12. Oscillations in Finite Quantum Systems, by G. F. Bertsch 
and R. A. Broglia, Cambridge University Press, Cam- 
bridge, England (1994). 

13. Optical properties of Metal Clusters, by U. Kreibig and M. 
Vollmer, Springer, Berlin (1995). 

14. H. Hovel, S. Fritz, A. Hilger, U. Kreibig, and M. Vollmer, 
Phys. Rev. B 48, 18178 (1993). 

15. B. N. J. Persson, Surf. Sci. 281, 153 (1993). 

16. Chaos in Classical and Quantum Mechanics by M.C. 
Gutwiller, Springer- Verlag, Berlin (1990). 

17. A. Kawabata, R. Kubo J. Phys. Soc. Jpn. 21, 1765 (1966). 

18. G. Bertsch, Comp. Phys. Comm. 60, 247 (1990). 

19. A. Rubio and L. Serra, Phys. Rev. B 48, 9028 (1993). 



1. J.-Y. Bigot, J.-C. Merle, O. Cregut, A. Daunois, Phys. 
Rev. Lett. 75, 4702 (1995). 

2. V. Halte, J. Guille, J.-C. Merle, J. I. Perakis, J.-Y. Bigot, 
Phys. Rev. B 60, 11738 (1999); J.-Y. Bigot, V. Halte J.-C. 
Merle, A. Daunois, Chem. Phys. 251, 181 (2000). 



